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Abstract 

The one-loop corrections to K — > 7r and K — > 2tt matrix elements of the 
electroweak penguin operator are calculated. General next-to-leading 
order relations between the K — > 7r and K —> 2n amplitudes are ob- 
tained. The fractional shift A 2 = 0.27 ± 0.27 is found for the 0(p 2 ) 
corrections to a recent chiral determination of {(tttt) j = 2 1 Q7.8 1-^°) ■ We 
explain why the sign for A2 is opposite to that expected from uni- 
tarization approaches based on the Omnes equation. We perform a 
background- field, heat-kernel determination of the divergent one-loop 
amplitudes for a general class of (V-A) x (V+A) operators. 



1 Introduction 



In the Standard Model, K — » 7T7r amplitudes are conveniently expressed in 
terms of the effective nonleptonic AS = 1 hamiltonian, 

Q 10 

(™\H A s=i\K) = -^=v ud v: s J2 <H(n) <HQ^>m , (i) 

where {Q«} are local four-quark operators, {c%(fJi)} are constants, and [i is a 
renormalization scale. In order to determine the K — * irir amplitudes, one 
must be able to calculate the matrix elements (7nr\Qi\K) ^. The study of 
such matrix elements has continued to be an active research area because of 
what it can teach us about the inner workings of low energy QCD. Interest 
has been recently heightened by the KTeV and NA48 announcements |l]] 
that e'/e ~ 20 -10 -4 , to be compared to the approximate theoretical relation 



- = 10 x 10" 4 



3.1 GeV- 3 - ((tttt) i=0 \Q 6 \K°) 



0.51 GeV" 3 - ((7r7r) /=2 |Q 8 |^ ) 



(2) 



evaluated in the NDR scheme at scale /i = 2 GeV. 

In the application of chiral symmetry to this problem, the K — > -kit 
matrix elements are expanded in powers of the external momenta and quark 
masses. It is especially attractive to consider matrix elements which do not 
vanish in the chiral limit, such as those of the electroweak penguin operators 

3 

Qi = 2 & di) v _ A Qq (?7 9i)v+A , 

~ q=u,d,s 

3 

Qs = ^(sidj) Y _ A }J Q q (qj qi) v+A , (3) 



E 

q=u,d,i 



where Q q is the electric charge of quark q and i, j are color labels. Working 
in the chiral limit, Donoghue and Golowich p| recently evaluated the leading 
chiral component to {(7T7t)/ = 2 | Qs\K°) and also ((7T7r)/ = 2 1 Q7I-K" ) . In this pa- 
per, we use chiral perturbation theory (ChPT) to extend this by calculating 
the 0(p 2 ) chiral corrections to matrix elements of Q7 and Qs- 

Given the flavor and chiral structure of Qj%, there exists a unique oper- 
ator at O(p ) that represents them in an effective chiral description, 

O ewp = gTr \X 6 UQU^} . (4) 
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The operator O ewp in turn belongs to a family of chiral operators O a b which 
transform as members of 8^ x 8r under chiral rotations, 



O ab = Tr 



(5) 



In Eqs. d),(D, A a is a Gell Mann matrix, Q = diag (2/3,-1/3,-1/3) is 
the quark charge matrix and U is the matrix of light pseudoscalar fields, 

U = exp(i\ k $ k /F) (k = l,...,8) , (6) 

where F is the pseudoscalar meson decay constant in lowest order. The 
coupling constant g in Eq. (|j) will depend on the 'parent' operator (Q7.8 
in our case) and can be obtained by comparison with an evaluation of the 
matrix elements performed in the chiral limit (analytically as in Ref. Q or 
QCD lattice-theoretic as in Ref. [||]). 

The standard first step in a ChPT analysis is to work in the chiral 
world of zero momentum and vanishing light quark mass. This leading 
chiral component, if nonzero, often provides the largest contribution to the 
matrix element. In the chiral limit, moreover, simple linear relations exist 
among amplitudes with differing numbers of Goldstone modes. In particular, 
knowledge of the K — > n matrix elements is sufficient to extract the physical 
K — » 7T7T matrix elements. 

Experience has shown, however, the necessity for calculating chiral cor- 
rections. This is especially true for amplitudes involving kaons, as the possi- 
bility of sizeable chiral corrections (of order 25% or more) cannot be excluded 
when passing from mx = to rtix = 0.495 GeV. @ This calls for an 0(p 2 ) 
analysis in the chiral representation of weak operators. At this order the 
chiral determination of K — > rrn matrix elements of Qr t s includes two basic 
ingredients: 

1. One-loop diagrams with insertion of vertices from O ewp (c/Eq. (Q)). 
We calculate in dimensional regularization, using a scale parameter 

2. Tree- level diagrams with insertion of the 0(p 2 ) counterterm operators 
(to be described in Sect. 2). 

In the K — > it sector, for example, the 0(p 2 ) amplitudes arise from the 
contributions depicted in Figure In Fig. 0(a), a loop begins and ends 

1 In Figs. strong and weak vertices are denoted respectively by large bold squares 
and small bold circles. 
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Figure 1: K — ► tt contributions: (a) weak vertex, (b) strong and weak 
vertices, (c) counterterm. 



at the electroweak vertex, whereas in Fig. |T](b) the electroweak vertex is an 
insertion in the loop which is produced by a strong interaction vertex. 

An important component of our calculation will be a leading-log estimate 
of the 0{p 2 ) corrections (c/Eq. (]4l|)). This provides only a partial estimate 
of the full 0(p 2 ) effect because it leaves out the contribution from the finite 
counterterms. Recall that divergences encountered in the loop diagrams 
are cancelled by divergences in the counterterm coefhcents, leaving just the 
finite loop and counterterm amplitudes. Each of these is dependent on the 
scale fx x , but the total amplitude is scale-independent. Unfortunately, the 
finite counterterm coefficients are unknown at this time. Therefore we base 
our estimation of 0(p 2 ) corrections on the finite part of the loop amplitude 
and determine the accompanying uncertainty by numerically studying its 
dependence on fj, x . 

In practice, first-principle calculations tend to reliably handle only the 
relatively simple K — > tt matrix elements. It is not clear that this infor- 
mation is sufficient to recover matrix elements for multipion final states. 
Indeed, reconstruction of the physical K — > tttt matrix elements is known to 
be nontrivial at next-to- leading order .0 An advantage of the ChPT approach 
is that it provides independent calculations of the K — > tt and K — > tttt tran- 
sitions. A significant aspect of our investigation will be to check whether 
the the 0(p 2 ) K — > tt matrix elements are sufficient to fully determine the 
the 0(p 2 ) K — > tttt matrix elements. 

Our presentation is arranged as follows. In Section 2 we list the 0(p 2 ) 
chiral operators representing Qj t s an d then perform the one-loop renormal- 
izaztion of the generating functional with insertions of the generic operator 
O a t, (see Eq. (||)). In Section 3 we present the 0(p 2 ) analysis for (tt\ Qt^\K) 

2 See Refs. 13, ]a, m for the analogous treatment of the 8_l x Ir and 27l X li? operators. 



3 



and {7nr\Qf t s\K) , including loop and counterterm contributions. Then, in 
Section 4 we perform a first numerical evaluation (based on chiral loops only) 
of the corrections to the leading chiral component of {tt^Oj ^\K) . The main 
results of this work are summarized and future studies are outlined in Sec- 
tion 5. Finally, we note that throughout the paper K — > nir matrix elements 
of a local operator O are expressed as 



M 



K- 



(7r( Pn )\0\K( Pk )) and M 



K — >7Tl7T2 



{ir( Pl )Tr(p 2 )\0\K( Pk )) 



(7) 



2 Counterterm Lagrangians 

In the following we first display the list of 0(p 2 ) counterterms which share 
the chiral properties of the O(p ) operator O ewp . Presentation of the various 
K — > -/r and K — ► tttt counterterm matrix elements is deferred to later in the 
paper. 

With each counterterm operator will be associated a coefficient which is 
a priori arbitrary. Some of these counterterm coefficients need to be singular 
in order to cancel divergences which appear in the one-loop amplitudes. In 
the second part of this section, we employ background field and heat kernel 
methods to determine the set of singular 0(p 2 ) coefficients. 

2.1 The Set of Counterterm Lagrangians 

There are seven effective operators at chiral order p 2 associated with the 
operator O ab of Eq. (|5|), 



with 



O 



O. 
O 
O 
O 



ab) 

L 

(ab) 



(ab) 



(ab) 
4 

(ab) 



Tr 
Tr 
Tr 
Tr 
i Tr 



ab 
ct 



(ab) 
k 



(8) 



k=l 



X a D^U\ b D^ 
KUDytA ■ Tr \\ b U ] D^U 
X a UX b D^D^U U ] + XaUD^D^U X b W 
X a UX b SU^ + X a USX b U^ 
X a U[P,Xb]U^ 
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o, 
o 



(ab) 



(ab) 



Tr 
Tr 





■ Tr 






• Tr 





(9) 



where the quantities S and P occurring respectively in and are 

defined as 



S = U ] y + Y ] U 



and 



P = i(u ] 



'x-x 



(10) 



The possible two-derivative dependence D^D^U has been eliminated by 
using the equations of motion, 



fill 



The case of interest for this paper (involving the operator O ewp of Eq. (|j)) 
is recovered with the replacements A a — ► Xq and A& — ► Q. Calculation of 
the set of counterterm amplitudes is straightforward and is summarized in 
Sect. O. 



2.2 Divergences in the One-loop Effective Action 

The effective chiral lagrangian on which the one-loop analysis is based is 



C = + O. 



str ~r ^ewp j 



where O ewp is defined in Eq. (|4j) and * s the familiar A 5 = lagrangian, 

£® = ^Tr(D„UD*rt) + !$-Tr(xrt + Ux*) ■ (13) 

We shall calculate quantum corrections about a solution U of the classical 
theory. To accomplish this, we employ background field and heat kernel 
methods.^] Thus we write 

(14) 



(12) 



U = U e iA 



where A = A a A a represents the quantum fluctuations. We begin by ex- 
pressing the generating functional at one-loop order as 



i J d*x[C(U)+£ ct (U)) . n du ^ e i J d*x [C(U)-C(U)} ^ ^ 



3 A summary of these techniques appears in Appendix B of Ref. 
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where C is given in Eq. (12) and C c t is the counterterm lagrangian of Eq. (||) 
with A a — > Xq and A;, — > Q. Keeping just the part quadratic in A, we have 



C(U) - C(U) = -^-A a (d^ + a + w) ab A b + ... 



(16) 



where 



jab 

v 

■ M 
aft 



Tttft 



Tr 



(7 



10 



.aft 



F 2 



[X a ,X b ] (u^dJJ + iUh^U + ir^ 



Tr 



A 6 f7 ( X a QX b + X b QX° 



1 



[[A a ,A 6 ] + ,Q] + )t/ 



(17) 



In the above and are source functions coupled to chiral currents. The 
integrand occurring in Eq. ([15]) is gaussian and thus allows for direct evalu- 
ation of the path integral, 



'1— loop 



d 4 x 



C{U)+C ct {U)\+-tx\og(d ll d^ + a + w) , (18) 



where 'tr' indicates a sum over spacetime coordinates as well as flavor labels. 
By using the heat kernel expansion, we identify the divergences in i?i_i 00 p 

as 



Z 



(div) 
1— loop 



1 



2(4vr) d / 2 



(Fx r 



where e 



4 - d and Tf v 



Tr 



^r^ + i(<T + u;) ; 



(19) 



[c^,c^p. There is a great deal of content 
in Eq. (|19|), including the one-loop divergent contributions to the Gasser- 
Leutwyler 0(p 4 ) strong lagrangian. g However, we require only the new 
piece which arises from the interference between the quantities a and w in 
Eq. pp, 



r(div) 
'1-loop 



2(4vr) d / 2 



d A x T - Tr [a w 



X 

J2 



d 4 x Tr [a w] , (20) 



where A is the singular quantity 



A = 



1 



16vr 2 



^-^(log^-T + l) 



(21) 
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Expressing our result in the operator basis (for the case A a — > A6 and At 
Q) of Eq. (P), we find for the general case of Nf flavors, 



-^l-loop — JP2 1 



20 2 + ^ (03 + 4 ) + Oe + O7 



(22) 



For three flavors (Nf = 3) as in this paper, we conclude that the one-loop 
generating functional will be finite provided we use 



(r) . Ac? 
- c y 

with 



c i = cf ) + ^d i (i = l,...,7) , (23) 



(24) 



di = dfe = , c?3 = c?4 = 3/2 , 
ci 2 = 2 , de = d 7 = 1 . 

3 K -> 7T and K -> 2tt Amplitudes at C(j9 2 ) 

When evaluated in the chiral limit, O ewp has the (tree-level) K ^ ir matrix 
elements 

7W {0) 7W (0) -0 f25) 

and the K — > 2-7T matrix elements 

These amplitudes have a very simple form, and it is not surprising that, as 
stated in Sect. 1, knowledge of the K — ► tt amplitudes is sufficient to yield 
the K — > 27r amplitudes. We state again that the numerical value of the 
coupling constant g can be deduced by referring to the work in Refs. [Q, [J. 
However, this information is not needed in the present analysis. 

The amplitudes at order p 2 will contain both loop and counterterm con- 
tributions, 

M (2) = M (loo P ) + M (ct) ^ (27) 

which are discussed separately in the subsections to follow. 
3.1 The One-loop Amplitudes 

Before proceeding to a discussion of the one-loop amplitudes, we consider the 
following technical matter. In the K — > it sector, unless the weak operator 
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is allowed to carry off nonzero four-momentum, the condition of energy- 
momentum conservation would be valid for physical states only in the SU(3) 
world of degenerate pseudoscalar meson masses (p 2 = p\ = fh 2 ). Therefore 
we shall allow the weak operator to transfer a four-momentum q. For the 
K — » nil matrix elements we can set q = 0, but we must keep a nonzero value 
of q for the K — > ir amplitudes. In particular, four-momentum conservation 
implies q = px — p^, where Pk,tt are four- momenta of the external kaon and 
pion. 

We now turn to the K-to-n matrix elements, for which we provide com- 
plete analytic expressions. This allows us to keep track of the q 2 (or equiv- 
alently the px ■ Pn) dependence, (ll]] and we obtain 

f -A(m 2 K ) + I K07r o(q 2 )] + V K o^o 



M 



(loop) 



M 



K°-^ n o — 9 
(loop) _ 

A'+^7T+ _ ' 



V2F 4 
' 2 1 
F K F n V2F 2 



34 A(m\) + 31 A(ml) + 9 A{m 2 n ) 



+ V 



(28) 



where g is the coupling constant of Eq. (||). In the above expressions, the 
functions I a f, are given by 



IkW (Q 2 ) = (3PK ■ Pn ~ m 2 K )B(q 2 ) + A{m 2 K ) - 

q 2 

I K +*+ (q 2 ) = (2m 2 K + ml- 3p K ■ p n ) B(q 2 ) + A{m 2 K ) 

+ g2+3(m f- m -W ), 



R(q 2 



where 



R(g 2 



m 2 K )B(q 2 ) + A(m 2 K )-A(ml) 



(29) 
(30) 



and A(m 2 ), B(q 2 ) are the regularized one-loop integrals 



iA(m 2 



iB(q 2 



lim 

d=A 



lim 

d=4 



2im 2 X + 



d d l 



1 



(27r) d I 2 



2i\ + fi 



i-d 



d d l 



(27r) d (l + q) 2 -m 2 K I 2 



rat 



(31) 



The corresponding divergent contributions for the K 
expressed in terms of A are 



7r amplitudes, as 



V 



8 



v K+ _ 



_\_ 

~p4 



(2 P K 



■Pit 



\Qm 2 K 



10 ml) 



(32) 



We consider next the two-pion matrix elements. Since one is ultimately 
interested in the physical on-shell result, we set q = 0. This allows us to 
replace cumbersome analytic expressions by their numerical values. At one- 
loop level the physical matrix elements will have both real and imaginary 
parts (as dictated by unitarity). Starting with the real parts, we have 



Tie 



-iM 



(loop) 



V2 



f( fl « +6 « log rlv) +I) « 



where the T>i are given for the K — ► irir amplitudes by 



V K o. 



V K ,_ 



A\/2 (Yi 



3 (m\ - mi) . 



The dimensionless coefhcients ctj, bi are collected in Table 1. 

Table 1 



Mode 




L h 


K U -► TT+7T- 


-1.195 


-1.300 


K° - 7T 7T 


-0.654 


-0.512 



Finally, the imaginary parts are found to be 



Im 



-?A/f (loop) 



V2g 1 p 



rn 



K 



-iM 



(loop) 



F K F% F 2 16vr 2 
y/2 g 1 /3 2 



F^F 2 F 2 16vr 



where /? = (! — 4m 2 /m 2 ^) 1 / 2 . 



3.2 The Counterterm Amplitudes 

Consider first the K — ► 7r transitions, for which we find 
,(ct) _ 4 



(33) 



(34) 



(35) 



F 2 



1 /4 \ 2 

g (ci - c 3 ) • + I -c 4 + c 5 + 2c 6 I 



9 



+ (c 4 + c 5 + c 6 ) ml 



M 



(ct) 



V2 
F 2 



1 2 \ 2 2 

-Cl + C 2 + 3 C 3 J Ptt • £>if - -^CmiK 



and next the K — > 7T7t counterterm amplitudes, 



_ju( cl ) 



x/2 

F 3 



g(ci - c 3 ) - 2c 4 - 2c 5 - 4c 6 ) m\ 



+ ( ~^( c i ~ c 3) - 4c 4 - 4c 5 - 2c 6 ) 



(ct) 



V% / 2 2 

^3 (m* - 



1 2 

— Cl — Co C 3 

3 3 



(36) 



(37) 



Upon adopting the prescription given in Eqs. (|23|) , (|24]) , the order p 2 ampli 



(2) 

tudes M\ are rendered ultraviolet finite. 



3.3 Construction of K — > rnr from K — > tt at 0(p 2 ) 

Let us consider now the finite parts of the counterterm amplitudes. We shall 
use the above results to show that knowledge of the K — > tt counterterm 
amplitudes implies knowledge of the K — > tttt counterterm amplitudes. As 
a result, each K —* tttt matrix element becomes expressible in terms of the 
K — > 7T7T loop amplitude and contributions obtined from the K — > tt sector. 

Recall from Eq. ( p3| ) that we denote the finite parts of the counterterm 
coefficients collectively as {cj r ^}. As a first example, we show how to fix 



the finite counterterm cf' which appears as a contribution to Ai 

(r) 



(ct) 

in Eq. (|37|). From Eq. (p6|), we can in principle obtain c^ 1 by taking a 
partial derivative with respect to the kaon squared-mass of the counterterm 
amplitude -M^}^ o- However, one works operationally with the full 0(p 2 ) 

Thus we have 



(2) 

amplitude M. K a 



j0 and the loop amplitude A^^^Lo- 



» 



2V2dm 2 x K °~ 



2V2dm 2 K V K °- 

(ct) 

The entire finite counterterm dependence of M K o_ 
manner, 



M 



(loop) 



(38) 



n n is obtained in like 



c\ + 3c 2 



2c 



(r) 



3F 2 



5 



\/2 9(p,r • pk) 



(2) 
K0- 



(loop) 



(39) 



10 



(ct) 

This leaves the set of remaining contributions to M K o_ 
find 



_ , for which we 



-i 



,(r) 



(r) 



3F 2 



9 



^2 

4 



_p2 

4 



4 #(?% • PK 

d d 



(m 



dm 2 
F 2 



dm? K 



M 



(2) 

K+- 

(2) 
K+- 



+ 







d 



y/2dm 2 K 
d 



M 



(2) 
K°- 



dm 2 K 



+ 



dm 2 

f 2 d_ 

2V2dm 2 K 



M 



(2) 
K+ 



(m 



(2) 
K°- 



M (loop) 



M {loop) 



M 



(loop) 



(40) 



The collection of finite counterterms in Eqs. (38)- ([5]) is seen to recover the 
entire content of the K — > tttt counterterm amplitudes. The appearance of 
partial derivatives in these equations suggests plotting the dependence of 
the K — > 7T amplitudes for each of the kinematic variables m|, m 2 , p n • px 
and numerically extracting the slopes. 



4 A First Estimate of the 0(p 2 ) Chiral Corrections 

At this point we have sufficient information to perform a first estimate of 
corrections to the physical matrix elements away from the chiral limit. Our 
estimate will involve only the chiral logarithms and, as discussed earlier, is 
therefore by no means complete. In particular, the answer in this approxima- 
tion depends upon the dimensional regularization scale \x x . An additional 
limitation is given by the fact that the chiral-log approximation does not 
distinguish between Qy and Qg. In both cases it gives the same fractional 
correction to the results obtained in the chiral limit. We turn now to the 
calculation of the fractional corrections and then compare our findings with 
those expected from an Omnes type unitarization approach. 



4.1 The Fractional Shifts A and A 2 

In order to recover the quantities commonly used in the literature, we must 
take matrix elements between K and a final two-pion state of definite isospin. 
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Denoting M.j = { (7T7t) j | Q7,s | -f^) , we can write each isospin amplitude, in- 
cluding chiral loops, as 

Mi = Mf ■ (1 + Aj) , 

where Aij is evaluated in the chiral world. We determine the central value 
of the 0(p 2 ) correction Aj by averaging the finite one- loop amplitude over 
the range fi x : 0.5 — * 1 GeV. The error bars in this determination are 
estimated from the range of values observed while varying the chiral scale 
[i x as indicated above. This uncertainty reflects our present ignorance of 
the low energy constants. Our procedure results in the 0(p 2 ) corrections, 

A = 0.98 ± 0.55 , A 2 = 0.27 ±0.27 . (41) 

These values omit possible contributions, of unknown magnitude, from the 
finite counterterms. It is clear the counterterm contribution cannot be very 
much smaller than that of the loop, else its scale dependence would not 
undergo the necessary cancelation against the loop component. In principle 
it could be much larger, in which case corrections exceeding the above ranges 
would occur. Even from our leading-log estimate, however, we see that 
sizeable deviations from the chiral limit predictions cannot be excluded and 
that a more detailed calculation of such corrections is called for. 



4.2 Comparison with Unitarization Procedures 

The possible impact of final-state interactions (FSI) on weak transition am- 
plitudes has recently drawn attention in the literature, p^, 13] Using the 



work of Ref. [12[] as an example, the K — > (tttt)j amplitude induced by the 
usual nonleptonic weak interaction is written in the form 

Ms) = Af\s - mDe^TZM , (42) 
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where is the derivative of Ai at the point s = m 2 and IZi(s) is a dis- 
persive corection factor. The above form is an Omnes-type solution which 
includes the effect of final-state rescattering. In the approximation of in- 



cluding just the two-pion final states, Ref. [12] obtains for the / = 0, 2 
rescattering corrections 

K (m 2 K ) = 1.41 ±0.06 , TZ 2 (m 2 K ) = 0.92 ±0.02 . (43) 

This means that, within errors, the effect of final state interactions enhances 
the 1 = amplitude by 41% and suppresses the 1 = 2 amplitude by 8%. It is 
revealing to compare these FSI effects with our 0(p 2 ) fractional corrections 
Ao and A2 given in Eq. (|4l"|). Our finding that Ao is large and positive is in 
qualitative accord with the above value for TZQ^m 2 ^), but our positive value 
for A2 conflicts with the suppression implied by TZ2{m 2 K ). 

The resolution of this apparent paradox lies in a more careful comparison 
of the two approaches. In Fig. ||(a)-(d), we display the four distinct one-loop 
contributions to our K — ► tttt electroweak penguin matrix elements. The 
contribution of Fig. |2](c) (the 'back-facing swordfish') is the only one with 
a two-pion intermediate state and is the only contribution whose iteration 
would occur in an Omnes-type resummation. A calculation of the fractional 
shift induced by Fig. |2](c) yields 

Aq (c) =0.36 ±0.12 , Af c) = -0.15 ±0.05 . (44) 



Both these values are in reasonable agreement with those implied by Eq. p3|). 
It is the other contributions in Fig. || (which we emphasize are demanded 
by the stringent requirements of chiral symmetry) which are dominant and 
which determine the overall sign of A2. 



5 Conclusion 

The calculation of matrix elements of local operators which enter the AS = 1 
weak effective hamiltonian is a crucial step in the phenomenology of nonlep- 
tonic kaon decays. In this paper we have focussed on some general aspects 
of the 0(p 2 ) corrections to the electroweak penguin operators 67,8 • Our 
results can be summarized as follows: 

1. Direct calculational methods (such as lattice-QCD, chiral sum rules, 
large-iV c QCD) usually deal with the simpler K — > ir transitions and 
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rely on chiral symmetry relations to infer the physical K — > tttt ma- 
trix elements. Although trivial in leading order (which is O(p ) in 
this case), these relations become more involved at next-to- leading or- 
der in the chiral expansion. We have systematically worked out the 
form of the K —* tt and K — > tttt matrix elements at 0(p 2 ) in the 
framework of ChPT, including one-loop diagrams and local countert- 
erms (as required by power counting) . We find that knowledge of the 
K — » 7r matrix elements is generally enough to fully determine the low 
energy constants entering in the K — > tttt matrix elements. Thus all 
chiral studies of this system, analytic or lattice-QCD, can safely use 
the K — * n mode to infer properties of the K — * tttt system. 

2. We have used background-field and heat-kernel techniques to perform 
the renormalization of the one-loop effective action in the presence 
of an insertion of the generic 8^x8^ operators O a b of Eq. (gj). Al- 
though the results given in Eqs. (|22l),(|24T) pertain specifically to the 
electroweak penguin operator of Eq. (||), the procedure is easily gen- 
eralized to applications involving other left-right operators. The di- 
vergent counterterm coefficients predicted by the general heat-kernel 
approach are in complete agreement with the divergences explicitly 
found in our one-loop calculation and thus provide a check on the 
calculation. 

3. Another important aspect of our work is the 'chiral-loop' estimate of 
the 0{p 2 ) correction to {{kit) i\Qi ,%\K) (cf Eq. (glj)). O n a qualita- 
tive level, we have answered the question of whether the corrections 
might be anomalously small (0(m 2 )) due to the absence of 0(m 2 K ) 
contributions. For example, recall that electromagnetic corrections to 
the K + — * 7r + 7r° matrix element of the octet weak lagrangian con- 
tains only small 0(am 2 T ) corrections even though the larger 0{am 2 K ) 
terms are allowed by power counting and symmetry arguments. We 
have found here that 0(m 2 K ) effects are not absent. This is why the 
leading-log result exhibits a strong dependence on the chiral scale [i x . 
Our analysis cannot, therefore, exclude the possibility of important 
corrections to the chiral determination. 

More quantitatively, we obtain the large 0(p 2 ) correction Ao = 0.98 db 
0.55 in the isospin-zero two-pion channel. It is not surprising that 
this agrees with Omnes-type resummation procedures [|^, [0| because 
strong rescattering effects are expected to play an important role for 
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the 1 = channel. Of greater interest is the 1 = 2 two-pion channel, 
since it bears on attempts to predict e'/e. Here, the situation is more 
subtle. The fractional correction A2 = 0.27 ± 0.27 is moderate in 
magnitude and positive in sign. The sign of A2 is opposite to that 
obtained in an Omnes-type resummation approach. As explained in 



Sect. L2 this is attributable to 0(p 2 ) contributions not included in the 
unitarized amplitudes. The underlying lesson is that unitarization of 
a part of the full amplitude can lead to valid results when rescattering 
is strong but can be misleading when rescattering is weak. 

This work represents a preliminary but necessary step towards a fully 
predictive analysis of the matrix elements (inr\Q7$\K) beyond the chiral 
limit. Work is underway to extend the dispersive analysis carried out in 
Ref. H to the general case of m quar k ^ 0. 
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